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INTRODUCTION 


In  the  last  five  years  the  increased  availability  of  hngh-speed 
electronic  computers  has  greatiy  enlarged  the  scope  of  data -reduction 
techniques  The  calculat'ons  of  double- layer  capac-ty  and  overpotential 
from  data  taken  using  a  dropping  electrode  in  conjunction  with  an  oscillo 
scope  are  not  difficult-  but  are  extremelv  red’cus  to  carry  out  with  an 
ordinary  desk  calculator  As  a  result  many  workers  have  not  obtained 
all  the  available  information  from  their  dam  and  have  relied  on  averaged 
values  which  often  vary  considerably  from  the  true  values 

This  report  describes  the  compu'er  programs  which  we  have 
developed  for  carrying  out  calculations  of  double  layer  capacity  and  over 
voltage  from  dam  taken  using  a  dropping  amalgam  electrode 

The  most  conveniently  available  computers  were  the  IBM  7094-1401 
comb.nation  which  was  used  for  long  calculations  ‘Program  3)  and  the 
IBM  1620  which  was  used  for  short  calculations  tFrcgrams  4  and  5)  The 
programs  were  prepared  m  the  FORTRAN  programming  language  since 
boi;h  of  these  computers  use  FORTRAN  compilers  as  standard  procedure 
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PROGRAM  3:  CAPACITY  OF  A  DROPPING  ELECTRODE  * 


The  major  program  discussed  here  is  one  written  for  the  IBM  7094, 
which  accepts  raw  capacity  data  and  calculates  the  double- layer  capacity, 
as  well  as  all  the  thermodynamic  parameters  which  can  be  obtained  from 
a  single  capacity- potential  curve.  Alternative  paths  are  provided  so  that 
the  capacity  can  be  calculated  without  performing  the  integrations,  or  so 
that  the  integrations  can  be  performed  on  capacity- potential  data  collected 
from  other  sources.  A  logical  diagram  of  the  program  is  shown  in  Fig.  1. 

The  experimental  technique  and  theory  on  which  the  calculations 
are  based  were  developed  primarily  by  Grahame  but  have  been  used 

by  many  other  workers. 

The  first  step  in  all  the  calculations  relating  to  dropping  electrode 
is  to  obtain  the  area  of  the  drop  at  the  moment  when  the  measurement  was 
taken.  A  convenient  form  for  the  dependence  of  area  on  time,  which  in¬ 
cludes  a  correction  for  the  back-pressure  due  to  interfacial  tension,  was 
proposed  by  Grahame  v  '  .  The  pressure  P  at  the  capillary  tip  results 
from  both  the  height  of  the  amalgam  above  the  tip  and  the  back-pressure 
due  to  interfacial  tension: 

P  =  ?gh  -  2(  Y/r)  (1) 

where  g  is  the  acceleration  of  gravity  (980.  6  dynes/g),  h  is  the  height 
of  the  column  of  amalgam,  p  is  the  density  of  the  amalgam,  V  is  its 
interfacial  tension,  and  r  is  the  radius  of  the  drop  at  the  end  of  the 
capillary. 

The  volume  flow  rate  is  taken  to  be  proportional  to  P  ,  and  the 
area  of  the  drop  is  obtained  by  integrating  the  volume  flow  rate  with 
respect  to  time,  assuming  a  spherical  shape  for  the  drop. 

The  assumption  of  a  spherical  drop  is  not  as  drastic  as  it  might 
appear,  since  as  long  as  the  drop  retains  its  cylindrical  symmetry,  the 
area/volume  ratio  is  independent  of  its  size.  The  error  resulting  from 


*  Programs  L  and  2  were  related  to  a  different  problem.  See  Tech. 
Memo.  #11:  "Hydrogen  Evolution  at  Solid  Indium  Electrodes.  " 
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the  normal  asphericity  of  the  drop  just  before  it  falls  can  be  avoided  by 
discarding  values  obtained  very  close  to  the  end  of  drop  lifetime.  The 
principal  uncertainty  in  area  occurs  when  the  drop  is  initially  formed 
and  ’s  not  yet  a  complete  sphere  This  problem  can  be  minimized  ex¬ 
perimentally  by  using  fine- pointed  capillaries  instead  of  the  conventional 
blunt  poiarographic  capillaries  ^  . 

Making  the  assumption  of  a  spherical  drop,  we  obtain  the  follow¬ 
ing  relation  for  the  area  of  the  drop  as  a  function  of  time: 

A  =  (36  7f  )’‘/3  (k  .  gh;)2/3  f,  .  3^_  f  \  */3j  ^ 

t  f  gh  V  3  m  t  /  ) 

where  k  is  a  proportionality  constant,  and  m  is  the  mass  flow  of 

amalgam  through  the  capillary,  which  is  independent  of  potential  ^  . 

The  proportionality  constant  can  be  evaluated  if  one  knows  the  drop  life¬ 
time  t*  ,  since  the  total  volume  of  a  drop  is  simp’y  m  t* /  ^  .  Evaluating 
(2)  at  the  lifetime  of  the  drop,  and  calculating  the  total  volume  of  the  drop, 
we  obtain  the  following  expression  for  the  proportionality  constant 


k.  =  ■?-  -  fi  -  UL  /  4 fff  i  13 

p  g  h  ^  g  h  ^  3  m  t*  / 

This  expression  may  be  substituted  in  (2)  to  give  an  equation  for  the  area 
of  a  drop  as  a  function  of  time  Note  that  when  h  is  very  large,  (2)  and 
(3)  together  reduce  to  the  usual  formula 


A  =  (36  7t  )1/3  (m  :/  £  )2/3 


(4) 


The  true  area  is  smaller  than  the  area  calculated  from  equation  (4)  be¬ 
cause  of  the  back-pressure 

Since  it  is  sometimes  difficult  to  obtain  accurately  reproducible 
drop  lifetimes,  some  sort  of  averaging  procedure  for  drop  lifetime  is 
desirable.  The  drop  lifetime  used  in  equation  (3)  is  calculated  from  the 
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interfacial  tension  V  and  an  average  value  for  the  diameter  a  of  the 
capillary  obtained  from  a  number  of  measurements.  The  drop  lifetime 
is  related  to  the  interfacial  tension  and  capillary  diameter  by  the  usual 
formula  for  the  drop- weight  measurement  of  interfacial  tension: 

t*  =  (  tf  Va/mg)  (5) 

Although  more  subtle  corrections  are  known  for  this  equation,  they 
are  net  conveniently  expressed  in  analytic  ferm,  and  were  not  necessary 
for  the  present  calculations. 

From  a  number  of  measurements  of  drop  lifetime,  an  average 
value  is  obtained  for  the  capillary  diameter.  This  can  be  checked  against 
actual  measurements  of  the  capillary  diameter  with  a  traveling  micro¬ 
scope.  From  this  value,  using  equation  (5)  ,  t*  is  calculated  for  each 
point,  and  used  in  (3)  to  obtain  the  proportionality  constant. 

The  advantage  of  this  averaging  method  is  that  only  a  few  drop 
lifetimes  need  be  actually  measured.  The  capacity  bridge  can  be  balanced 
at  the  same  point  in  the  oscilloscope  trace,  whether  or  not  the  end  of  the 
drop  actually  appears  on  the  oscilloscope  screen  . 

In  the  above  equations  a  value  must  be  assumed  for  the  interfacial 
tension  in  the  initial  phase  of  the  calculation,  but  once  the  capacity  data 
has  been  integrated,  a  much  better  value  for  the  interfacial  tension  is 
known.  To  provide  the  most  rapid  convergence,  the  best  possible  starting 
values  for  interfacial  tension  are  required  These  are  obtained  in  the 
present  calculation  by  intergrating  an  assumed  quadratic  form  for  the 
integral  capacity,  which  may  be  obtained  from  literature  data,  oi  guessed 
by  analogy  with  similar  systems.  The  exact  values  used  affect  only  the 
rapidity  of  convergence,  not  the  final  answer. 

The  integral  capacity  is  assumed  to  be  given  by  the  quadiatic  form 

K  =  Ko  +  K,  +  K2  ^  2  (6) 


-  4  - 


where  V  is  the  rational  potential  (potential  with  respect  to  the  point  of 
zero  charge).  The  interfacial  tension  is  then  given  by 

Vmax-10{-^  ¥  2  +  ¥3  +  ^  ¥  4  ' 

z  o  4  J 

where  ^max  is  the  maximum  value  of  the  interfacial  tension,  which 
must  be  measured,  obtained  from  the  literature,  or  estimated  by  correla¬ 
tion  with  similar  systems.  For  mercury  in  dilute  aqueous  solutions, 

Y  max  =  426'2  ei'g  cm"2,  KQ  =  27.0  |dfcm'2,  Kr  =  18.65  )4f  cm"2 

volt  ,  and  K2  =  9.65  (4f  cm  2  volt  2  .  The  factor  10  is  required  111 

equation  (7)  because  of  the  conventional  units,  one  uf  cm'2  volt  2  is 
_  9 

equal  to  0.  1  erg  cm 

Further  details  regarding  the  relation  of  drop  area  to  drop  lifetime 
may  be  found  in  the  literature  . 

Once  the  area  of  the  drop  at  the  point  cf  balance  has  been  calculated, 
the  values  of  capacity  and  series  resistance  read  from  the  bridge  can  be 
converted  to  the  desired  quantities:  capacity  per  unit  area  and  polariza¬ 
tion  resistance.  The  capacity  per  unit  area  is  simply 


C 


Cexp  /  A 


(8) 


The  polarization  resistance  is  obtained  after  ccrrecting  the  total  measured 
res. stance  for  the  resistance  of  the  electrolyte  and  the  resistance  of  the 
amalgam  in  the  capillary.  If  a  cylindrical  or  spherical  counter  electrode 
of  radius  ^  surrojnds  the  drop,  then  the  resistance  R  cf  the  solution  is 
given  by 


R  =  — 1 1 _ 

KA(r+)  ) 


(9) 


where  K  is  the  specific  conductance  of  the  electrolyte,  r  is  the  radius 


arid  A  is  the  area  cf  the  drop  at  the  balance  point.  The  polarization 
resistance  is  given  by 


A 


R 


R 


cap 


(10) 


where  R  is  the  resistance  measured  on  the  bridge,  R  is  given  by 

6Xp 

equation  (9),  and  Rcap  is  the  resistance  of  the  thread  of  amalgam  in  the 
capillary.  Rcap  may  be  calculated  from  the  dimensions  of  the  capillary 
and  the  resistivity  of  the  amalgam.  It  is  rarely  more  than  one  or  two  ohms. 

The  next  portion  of  the  calculation  involves  integration  of  the 
measured  differential  capacity  values  to  obtain  the  surface  charge,  integral 
capacity,  and  interfacial  tension. 

Measurements  of  capacity  are  usually  taken  at  quite  closely-spaced 
intervals  of  potential  (0.01  to  0.05  volts).  Since  the  values  may  scatter 
slightly  about  a  smooth  curve,  a  simple  trapezoidal  rule  was  used  for  the 
integration  instead  of  a  more  sophisticated  numerical  integration  algorithm. 
The  surface  charge  is  obtained  by  integrating  the  differential  capacity: 


q  = 


C  dif 


(10 


where  (f  is  the  rational  potential.  A  second  integration  gives  the  inter- 
facial  tension 


V  =  ^  max  '  10  J  <**1 


(12) 


The  absolute  value  of  V  depends,  of  course,  cn  the  value  of  Y  max  , 
which  must  be  obtained  independently.  The  integral  capacity  is  obtained 
from  the  surface  charge: 

K  =  q/  Vjf  U3) 
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V 


Once  these  values  have  been  calculated,  the  interfacial  tension  is 
used  in  the  first  part  of  the  program  to  recalculate  the  area  of  the  drop, 
and  better  values  of  the  differential  capacity,  surface  charge,  integral 
capacity,  and  interfacial  tension  are  obtained.  The  iteration  is  continued 
until  successive  values  of  interfacial  tension  agree  within  0.01  erg/cm^  , 
The  final  values  do  not  depend  on  the  initial  assumption  for  the  capacity 
function  (eq.  6);  this  affects  only  the  rapidity  of  convergence. 

The  next  step  in  the  calculation  is  to  find  the  potential  of  the  outer 
Helmholtz  plane  and  to  calculate  the  contribution  to  the  capacities  from 
the  diffuse  double  layer.  The  potential  of  the  outer  Helmholtz  plane,  'f', 
is  given  by  the  equation 


K°  (  (f  -  ^  ) 


2  D  Dq  RT  C 
7t 


sinh 


£_! 

2  RT 


(14) 


where  K°  is  the  integral  capacity  of  the  Helmholtz  double  layer,  D  is 

the  dielectric  constant  of  the  solvent,  DQ  (1.  1128  x  10  coul  volt  cm  ) 

is  a  dimensional  constant,  R  is  the  gas  constant,  C  is  the  concentration 

of  a  1-1  electrolyte,  T  is  the  absolute  temperature,  and  F  is  the  Faraday 

constant.  Analogous  equations  have  been  derived  for  electrolytes  other  than 

simple  univalent  ones,  ^),  and  treatments  have  been  extended  beyond  the 
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simple  Debye-Huckel  theory  v  '  . 

With  the  constants  evaluated  and  combined,  eq.  (14)  becomes 


K°  (  y  -  ^  )  =  11.72 


This  equation  is  most  easily  solved  by  iteration  using  Newton's  approxima¬ 
tion  method. 

Initially,  since  we  do  not  know  K°  ,  we  approximate  K°  by  K  , 
the  total  integral  capacity,  and  calculate  a  provisional  value  for  "vf’  , 
using  equation  (15)  .  Then,  using  this  approximate  value  for  "'f'  ,  the 
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integral  capacity  of  the  diffuse  double  layer,  Kd  ,  is  calculated: 

Kd  ■  ¥*  fiHr  “  (“  “  H?)} 

The  integral  capacity  of  the  Helmholtz  layer  is  obtained  from  the  total 
integral  capacity  and  the  capacity  of  the  diffuse  double  layer: 

_L  =  J_  -  J_ 

K°  K  Kd 


(16) 


(17) 


The  improved  value  of  Kq  is  then  used  in  (15)  to  obtain  a  better  approxima¬ 
tion  to  ‘vf>  .  The  iteration  is  continued  until  two  successive  values  of  Kd 

9 

agree  to  within  0.  01  jmf/cm  . 

The  differential  capacity  of  the  diffuse  double  layer  is  calculated 
from  the  equation: 

Cd  =  228.5  (^*\  J  C  — '  cosh 
\  T  /  T  298 

and  the  differential  capacity  of  the  Helmholtz  double  layer  obtained  from 
this: 


J_  =  J_  -  1 
c°  c  c( 


(19) 


This  iteration  often  converges  after  less  than  five  cycles.  The 
convergence  is  most  rapid  for  concentrated  solutions,  since  then  Kd  is 
very  close  to  K  . 

In  Appendix  I  is  given  a  complete  listing  of  Program  3,  together 
with  the  required  forms  for  the  input  data,  and  sample  input  and  output 
data. 
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PROGRAM  4-  OVERVOLTAGE  MEASUREMENTS  AT  A  DROPPING 
_ _  ELECTRODE 


Conventionally,  when  cur  rent -voltage  curves  are  measured  using 
a  dropping  electrode,  and  a  simple  potentiometer  polarizing  circuit, 
either  the  time-average  current  or  the  current  at  the  end  of  a  drop  life¬ 
time  is  used  to  calculate  the  current  density.  Such  a  method  is  adequate 
if  the  currents  are  very  small,  or  if  only  approximate  current-potential 
curves  are  required.  However,  under  conditions  where  currents  of  the 
order  of  1  milliampere  are  to  he  drawn,  or  where  more  accurate  measure¬ 
ments  of  overvoltage  are  required,  a  better  method  is  to  observe  the 
instantaneous  value  of  the  dropping  electrode  potential,  while  applying  a 
constant  current  to  the  cell  through  an  external  polarizing  circuit.  The 
difference  between  the  drop  potential  ar.d  a  known  bias  potential  can  be 
displayed  on  an  oscilloscope,  and  under  favorable  conditions,  the  instan¬ 
taneous  potential  of  the  drep  can  be  measured  to  within  1  or  2  mv.  A 


typical  trace  :s  shown  in  Pig.  2.  The  constant  current  ran  lx?  read  from 
an  accurate  microamrneter,  ar.d  the  current  density  calculated  once  the 
drop  area  is  known. 

a  he  calculation  of  drop  area  is  essentially  the  same  as  we  have 
described  for  the  capacity  measuiements-  .however,  for  overvoltage 
measurements  the  accuracy  required  ’s  net  so  great.  An  error  of  1  mv 
m  measuring  hydrogen  overvoltage  is  equivalent  to  an  error  of  2%  in 


cu  rent  density'.  This  means  that  the  back-pressure  correction,  which  n 
significant  for  capacity  measurements,  i?  unimportant  for  overvoltage 
measurements  and  can  lx?  dispersed  with.  Thus  the  stmpiem  equation  (4) 
can  be  used  for  calculating  the  area. 


A  (367 T  }i/3  (m  t,  {  )2/3 


(4) 


lo  ease  the  tedium  of  measuring  several  potential- time  points  on 
each  of  several  hundred  oscilloscope  photographs  to  within  0.  01  inches, 
a  Gerber  Analog-ro-lLgital  converter,  or  graph  reader, was  used  This 
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device  reads  the  x  and  y  coordinates  of  a  manually-set  cross-hair, 
and  enters  the  information  on  a  punched  card  in  a  pre -determined  format. 
Because  the  coordinates  are  given  on  an  internal  integral  scale.  200  units 
per  inch,  it  was  necessary  to  include  provision  for  horizontal  and  vertical 
calibrations  in  the  program. 

A  complete  listing  of  Program  4,  together  with  the  forms  for  the 
input  data,  and  sample  input  and  output  data,  are  given  in  Appendix  II. 
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PROGRAM  5:  LEAST-SQUARES  ANALYSIS  OF  CURRENT-POTENTIAL 
_ CURVES  _ _ 

The  punched  card  output  from  Program  4  was  used  in  another  program 
which  fitted  current-density  i  and  overpotential  7|  measurements  to  a 
theoretical  curve  derived  by  Frumkin  ^  on  the  assumption  that  dis¬ 
charge  of  an  ion  is  the  slow  step  in  the  electrode  reaction  mechanism. 


7  ■ 


(X  F 


in  (  i  /  i0'  ) 


(20) 


where  V'  is  the  potential  of  the  outer  Helmholtz  plane,  calculated  according 
to  equation  (14).  The  integral  capacity  is  assumed  to  be  given  by  a  quadratic 
expression  (eq.  6).  Eq.  20  can  cast  into  the  form  of  a  straight  line 


y 

=  0(  x  +  B 

(21) 

where 

X 

(22) 

y 

■<('+—  L  i 

F 

(23) 

B 

^  ~  h  i0 

F  0 

(24) 

The  values 

of  x  and  y  are  calculated, 

and  a  least-squares  fit  made  to  a 

straight  line.  The  best  values  of  the  parameters  o(  and  -  log  i(  are 
calculated,  together  with  their  95%  confidence  limits. 

To  fit  the  current  density-potential  data  to  an  ordinary  Tafel  equation 

1  =  ~r  <  i/i0 )  (25) 


it  is  necessary  only  that  'f  should  be  zero;  this  is  accomplished  by  setting 
all  the  coefficients  of  the  quadratic  expression  for  capacity  (eq  6)  equal  to 
zero 
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A  complete  listing  of  Program  5,  together  with  the  requisite  forms 
for  the  input  data,  and  sample  input  and  output  data,  is  given  in  Appendix  III. 

CONVERSION  OF  1620  PROGRAMS  FOR  7090 

Conversion  of  a  program  written  for  the  IBM  1620  for  use  on  the 
IBM  7090  or  7094  involves  only  the  modification  of  the  input-output  state¬ 
ments.  To  convert  Programs  4  and  5  for  use  on  the  7090,  change  every 
"READ"  statement  to  "READ  INPUT  TAPE  5"  ,  and  every  "PUNCH"  state¬ 
ment  to  'WRITE  OUTPUT  TAPE  6".  Tiic  only  exception  is  statement  53 
of  Program  4,  which  must  produce  punched  card  output,  and  which  should 
read  "WRITE  OUTPUT  TAFE  9". 

Under  some  monitor  systems,  all  of  these  changes  may  not  be 
necessary. 

CONVERSION  OF  7090  PROGRAMS  FOR  1620 

Conversion  of  a  7090  program  for  use  on  the  1620  is  often  more 
difficult  than  the  reverse,  and  indeed  may  be  impossible  because  of  the 
more  limited  storage  facilities  of  the  smaller  computer.  The  7090  program 
described  here  (Program  3)  requites  no  more  than  4000  words  of  memory, 
and  so  can  lie  adapted  to  the  1620  with  only  slight  modifications.  These 
modifications  are  as  follows 

1.  Change  all  statements  'READ  INPUT  TAPE  5"  to  "READ". 

2.  Change  all  statements  "WRITE  OUTPUT  TAPE  6"  to  "PUNCH". 

3.  Change  all  statements  "WRITE  OUTPUT  TAFE  9"  to  "PUNCH", 
(the  punched  output  for  re-entry  to  the  program  will  have  to  be 
separated  from  that  which  is  to  be  printed. ) 

4.  Instead  of  the  REREAD  subroutine,  the  title  card  must  be  read 
as  a  dimensioned  variable.  Delete  the  REREAD  statement  at 
the  beginning  of  the  program.  Include  in  the  dimension  state¬ 
ments,  a  variable  TITLE  (20). 
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Statements  1  and  2  should  read: 

1  READ  2,  KK,  TITLE 

2  FORMAT  (II,  13,  1914) 

Statements  350  and  351  should  be  deleted. 

5.  The  functions  defined  in  statements  123,  38,  39,  179,  and 
180  at  the  beginning  of  the  program  must  be  included  ex¬ 
plicitly  in  the  statements  which  refer  to  them:  statements 
28,  50,  152,  184,  185,  191,  197,  198,  205,  211. 

Unless  only  a  very  small  amount  of  data  is  to  be  processed,  or  a 
faster  computer  is  not  available,  it  is  not  recommended  that  Program  3 
be  used  on  the  1620,  since  the  time  for  compiling  and  executing  it  will 
be  inordinately  long. 
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APPENDIX  I 


Figure  1  -  Logical  diagram  for  Program  3 
Listing  of  Program  3 
Text  describing  input  data  for  Program  3 
Sample  data 
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LOGICAL  DIAGRAM  FOR  PROGRAM  3  (CAPACITY) 


(  OUTPUT  DATA  ^ 


1 


CAPACITY  OF  A  DROPPING  AMALGAM  ELECTRODE  -  PROGRAM  3 


5/09/64 


C2 

C3 

C4 

C 

E 


123 

38 

39 

179 

180 
C 

1 

350 

351 
2 

313 
420 

314 
115 
316 


INPUT  OFF-LINE  ON  A2  -  LOGICAL  5  3 
PRINTEO  OUTPUT  OFF-LINE  ON  A3  -  LOGICAL  6  3 
PUNCHED  OUTPUT  OFF-LINE  ON  A5  -  LOGICAL  9  3 

REREAD  3 
DIMENSION  NUMK80),  NUM2I80),  CAPG  (80),  ATENS  (80)  3 
DIMENSION  ETAK80),  Tl(BO),  ETA(RO),  TIME(80),  CAPX(80)»  RESX  (80)  3 
DIMENSION  AKEA( 80 ) »  CAPI80),  RES(80),  PHI(80)»  CHG ( 80) , CAPI NT ( 80 )  3 
DIMENSION  CAPH( 80 ) *  CAPIGI80),  CAPIH(80),  TENSIBO),  PSI (80)  3 
GAMMAF  (P)  *  P»P*(CAP0/2.  ♦  CAPi*P/3.  ♦  CAP2*P»P/4.)  *  10.  3 
8RAFI  T»G)  *  (  l.-(Q*G)/( T**( 1./3. ) ) )**(2  ./3. )  3 
TSTAKF(G)  *  (UtAM»G)/(AM*0. 31214)  3 
S INHF ( / )  »  0.5(EXPF(Z)-EXPF(-ZI)  3 
COSHF(Z)  *  0.5(EXPF(Z)*EXPF(-Z))  3 

READ  INPUT  TAPE  5,  350,  KK  3 
FORMAT  (III  3 
READ  INPUT  TAPE  99,  2  3 

FORMAT  ( 72H - THIS  SPACE  IS  FOR  A  TITLE - 3 

1 - ,  3 

IF  (KK  -  1)  314,  316,  420  3 
IF  (KK  -  2)  314,  316,  314  3 
WRITE  UUTPUT  TAPE  6,  2  3 
GO  TO  l  3 
WRITE  OUTPUT  TAPE  6,  2  3 


124  REAC  INPUT  TAPE  5,5,RHU,  HGT ,  AM,  DELTA,  TEMP  3 

5  F0RMAT(4F8.3,F8.2)  3 

255  AMT  »  AM  3 

230  IF  (AMT  -  999.99)  6,  231,  232  3 

232  WRITE  OUTPUT  TAPF  6,  233,  AMT  3 

233  FORMAT  (/51H  THE  TEST  NUMBER  FOR  MASS  FLOW  RATE  IS  WRONG.  AMT  =  3 

1  F8.3)  3 

324  GO  TO  1  3 

231  READ  INPUT  TAPE  5,  234,  AMI,  AM2  3 

234  FORMAT  (2F8.3)  3 


2 

3 

4 

4A 

5 

6 

7 

8 

8  A 
88 
8C 
8D 
8E 

9 

9  A 
9B 

10 
11 

1 1  A 
12 
12A 
12B 
12C 

13 

14 
14A 
148 
14C 
14D 
14  D 
14E 
14F 
14G 


C 


C 


C 


6 

READ  INPUT  TAPE  5,7,  GMAX,  ETAZ, 

CONC , SKA  PPA , 

CRCS 

3 

15 

7 

FORMAT  (F8.2,  FB.3,  2F8.5,  F8.3) 

3 

16 

300 

IF  (RHO  -  99.999)  8,  19,  301 

3 

17 

301 

WRITE  OUTPUT  TAPE  6,  302,  RHU 

3 

17  A 

302 

FORMAT  I/51H  THE  TEST-VALUE  FOR 

RHU  IS 

WRONG. 

IT  IS  NOT  99.999-  / 

3 

178 

115H  INSTEAD,  IT  IS  HO. 5) 

3 

17C 

320 

GO  TO  l 

3 

17D 

19 

N  *0 

3 

18 

20 

N  =  N  ♦  1 

3 

19 

121 

READ  INPUT  TAPE  5,  21,  ETA(N),  CAPIN), 

NUM1IN) 

,  NUM2IN) 

3 

20 

21 

FORMAT  (F8.3,  8X ,  F8.2,  BX,  214) 

3 

21 

122 

IF  (ETA(N)  -  3.333)  20,  26,  23 

3 

22 

23 

IF  (ETA(N)  -  6.666)  24,  26,  24 

3 

23 

24 

WRITE  OUTPUT  TAPE  6,  17,  ETA(N), 

CAPIN) 

3 

24 

321 

GO  TU  1 

3 

24  A 

26 

N  =  N-l 

3 

25 
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126  WRITE  OUTPUT  TAPE  6,  27,  TEMP,  CONC,  ETAZ,  GMAX  3 

27  FORMAT  (//8H  TEMP  *  F8.2/  8H  CONC  ■  F8.5/  8H  ETAZ  *  F8.3/  8H  GMAX  3 

I»  F8. 2  )  3 

127  GO  TO  70  3 

C 


8  REAO  INPUT  TAPE  5,4,  CAPO,  CAP1,  CAP2  3 

4  FORMAT  (3F8.3)  3 

9  WRITE  OUTPUT  TAPE  6,  10,  RHO,  GMAX,  CAPO,  HOT,  ETAZ,  CAP1,  AM,  3 

1C0NC,  CAP2,  DELTA, SKAPPA,  TEMP,  CRES  3 

10  FORMAT (//8H  RHO  ■  FR.3,6X,8H  GMAX  ■  F8.2,6X,8H  CAPO  ■  F8.3/8H  HGT  3 
1  *  F8.3,6X,8H  ETAZ  «  F8.3,6X,8H  CAP1  >  F8.3/8H  FLOW  >  F8.3,6X,8H  3 

2C0NC  «  F8.5,6X,8H  CAP2  «  F8.3/8H  DELTA-  F8.3,6X,10H  KAPPA  «  F8.6,  3 

34X,  8H  TEMP  >  F8.2  /  8H  CRfcS  >  FB.3)  3 

C 

250  IF  (AMT  -  999.99)  11,  251,  256  3 

256  WRITE  OUTPUT  TAPE  6,  233,  AMT  3 

257  GO  TO  1  3 

251  WRITE  OUTPUT  TAPE  6,  252,  AMI,  AM2  3 

252  FORMAT  (/58H  THE  MASS  FLOW  RATE  USED  IS  A  LINEAR  INTERPOLATION  BET  3 

1WEEN  /  7H  AMI  «  F8.3,  llH  ANO  AM2  >  F8.3)  3 

C 


C46  WE  NOW  EVALUATE  THE  CAPILLARY  01 AME TER 

11  M  »  0 

12  M  ■  M  ♦  1 

13  REAC  INPUT  TAPE  5,  14,  ETA1(M),  T1(M) 

14  FORMAT  (F8.3,  F8.5) 

15  IF  (ETAKM)  -  3.333)  12,  18,  16 

16  WRITE  OUTPUT  TAPE  6,  17,  ETAl(M),  T 1 ( M > 

17  FORMAT  ( / 58H  THERE  IS  NO  PROPER  ENO-OF-DATA  CARD  -  THE  LAST  CARD  R 
IE AD  F8.3,2X,F10.5> 

118  GO  TO  1 

18  M  ■  M  -  1 

119  WRITE  OUTPUT  TAPE  6,  22 

22  FORMAT  (//20H  OROP  LIFETIME  OATA  /32H  POT  TIME  GAMMA  Dl 
IA  /» 

32  FM  -  M 
139  NM  ■  l 
125  SOI  A  -  0 
25  00  31  I  »  1 , M 
303  IF  (NM  -  1)  28,  28,  330 

330  DO  332  II  *  1,  N 

331  IF  (ETA(II)  -  ETAKI))  332,  333,  333 

332  CONTINUE 

333  GAM  *  TENS!  II  ) 

305  GO  TO  236 

28  GAM  «  GMAX  -  GAMMAF( ETA l( I )-  ETAZ) 

400  ATENS(I)  ■  GAM 

236  IF  (AMT  -  999.99)  129,  237,  238 

238  WRITE  OUTPUT  TAPE  6,  233,  AMT 

239  GO  TO  l 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


237 

FI  » 

I 

240 

AM  * 

AMI 

♦ 

( AM2  -  AM 

1 ) 

*  FI/FM 

129 

01  A  = 

0. 

312 

14*  AM  • 

Tl( 

I)/GAM 

130 

SD I A 

-  SD  IA 

*  Dl  A 

306 

IF  (NM  - 

1) 

29,  29, 

31 

26 

27 

28 

29 

31 

30 

32 

33 

34 

35 

36 

37 

37A 

37B 

37C 

370 

37E 

37F 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 
53A 

53  B 
53C 

54 
54A 

54  B 
54C 
540 
54E 

54  F 
54G 

55 

55  A 
55B 
55C 

550 

55E 

56 

57 
57A 
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29  WRITE  OUTPUT  TAPE  6,  30,  ETAKI),  TI(I),  GAM,  OIA 

30  FORMAT  (F8.3,  F8.5,  F8.2,  F8.5) 

31  CONTINUE 

33  01  AM  ■  SOIA/FM 

307  IFINM  -  I)  34,  34,  242 

34  WRITE  OUTPUT  TAPE  6,  35,  OIAM,  M 

35  FORMAT  (/33H  AVERAGE  DIAMETER  OF  CAPILLARY  *  F8.5,I2H  MILLIMETERS/ 

126H  NUMBER  OF  MEASUREMENTS  *  12) 

C 

C48  THE  NEXT  THING  TO  00  IS  CALCULATE  THE  AREA  AND  THE  CAPACITY 

242  IF  (AMT  -  999.99)  36,  309,  244 
244  WRITE  OUTPUT  TAPE  6,  233*  AMT 
136  GO  TO  1 

36  COA  *  4.835974*1 ( AM*. 001/RH0)**! 2. /3. ) ) 

37  Q  *  0 .00493164/1 HGT * ( AM* .00l*RH0*RH0) •* ( 1 . / 3 • ) ) 

309  IF  (NM  -  l)  57,  57,  150 

57  WRITE  OUTPUT  TAPE  6,  2 
401  WRITE  OUTPUT  TAPE  6,  58 

58  FORMAT! //57H  FIRST  APPROXIMATION  BASED  ON  ASSUMED  INTERFACIAL  TENS 

1I0N  / /66H  RUN  POT  TIME  EXP-CAP  EXP-RES  AREA  C 

2AP  PUL-RES  /66H  NO.  VOLTS  SEC  M-FARADS  OHMS  SO  MM 
3  MF/SQCM  0HM-CM2/ ) 

C 

40  N  *  0 

41  N  *  N  ♦  1 

42  READ  INPUT  TAPE  5,  43,  ETA(N),  TIME(N),  CAPX(N),  RESX(N),  NUMl(N), 
1  NUM2IN) 

43  FORMAT  (F8.3,  F8.5,  2F8.3,  214) 

44  IF  (ETA(N)  -  3.333)  41,  47,  45 

45  IF  (ETA(N)  -  6.666)  46,  47,  46 

46  WRITE  OUTPUT  TAPE  6,17,  ETA(N),  TIME(N) 

322  GO  TO  1 

47  N  »  N-l 

246  FN  «  N 

48  00  50  I  «  1,N 

49  PHI ( I  )  »  ET All)  -  ETAZ 

50  TENS!  I)  *  GMAX  -  GAMMAF  (PHI(D) 

C 

150  TEST  *  TENS(l) 

51  DO  55  I  *  l»N 

243  IF (AMT  -  999.99)  152,  247,  152 

247  FI  *  I 

248  AM  «  AMI  ♦  ( AM2  -  AMI)  •  FI/FN 

253  COA  «  4.835974*1 ( AM* .OOl/RHO ) •*( 2./3. ) ) 

254  0  *  0. 00493164/1 HGT*( AM*. O01*RHO*RHO)**(l./3.)) 

152  RATIO  -  BKAF  ( T IME( I ) , TENS! I ) )/BRAF ( TSTARFI TENS ( I ) ) , TENS ( I ) ) 

153  ARE A( I )  «  COA  •  RATIO  •( T IMEI I )•• ( 2. /3 .  ) ) •  100. 

52  CAP ( I  )  *  CAPXI  D/IAREAI  I)*.01) 

53  RAO  «.0282095*SQRTF( AREA! I ) ) 

54  SOLR  *  DELT A*RAD/ (  SK APPA* ( RAD  *  DELTA)) 

55  RES ( I  )  *  O.Ol»AREA( I )•( RE SX( I)-CRES)  -  SULR 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


58 

59 

60 
62 

62  A 

63 

64 

65 

67 
67A 
67B 
67C 

68 
69 

72A 

72B 

73 

74 

75 

76 

77 

76 

79 

80 
80 
81 
82 

83 

84 
84A 

85 
85A 

86 

87 

88 

89 

90 
90A 
90  B 
90C 
90D 
90E 

91 

92 

93 

94 

95 

96 


THE  FIRST  TIME  AROUND  WE  PRINT  THE  RESULTS  AND  SEE  IF  WE 
ING  TO  INTEGRATE.  NEXT  TIME  AROUND  WE  SKIP  THAT  PART. 

56  IF  (NM  -  1)  163,  60,  59 

60  WRITE  OUTPUT  TAPE  6,  61 , ( NUM 1(1),  NUM2II),  ETA(l),  T I  ME ( I ) 


ARE  GO- 


, 


3 

3 

3 

3 


97 

98 

99 
100 
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; 

1  (I).  RESXU),  AREA!  I ) .  CAP(I),  RES(I),  1 

-  1,  N) 

3 

101 

61 

FORMAT  (1X.2I4.  2X»  F7.3,  F8.5,  5F8.3) 

3 

102 

159 

r 

GO  TO  63 

3 

103 

\ 0 

59 

IFIPHIIJ))  62.  73,  76 

3 

104 

62 

WRITE  OUTPUT  TAPE  6,  94,  PHI(J),  J 

3 

105 

162 

r 

GO  TO  1 

3 

106 

l# 

163 

WRITE  OUTPUT  TAPE  6,  164,  NM 

3 

107 

164 

FORMAT  I/41H  THE  ITERATION  INDEX  IS  LESS  THAN  1,  NM  •  13) 

3 

108 

165 

r 

GO  TO  1 

3 

109 

U 

63 

IF  IETAIN+1 )  -  3.333)  65,  70,  64 

3 

110 

64 

IF  I ETAI N+l )  -  6.666)  65,  67,  65 

3 

111 

65 

WRITE  OUTPUT  TAPE  6,  66,  ETAIN+1) 

3 

112 

66 

FORMAT  I/46H  WE  HAVE  LOST  THE  END-NUMBER . 

INSTEAD  WE  HAVE  FR.3) 

3 

113 

166 

GO  TO  1 

3 

114 

67 

WRITE  OUTPUT  TAPE  6,  68,  N 

3 

115 

68 

FORMAT  I/52H  THE  TEST  NUMBER  WAS  6.666  -  NO  INTEGRATION  WAS  DONE/ 

3 

116 

112H  THERE  WERE  12.  25H  DATA  POINTS  IN  THIS 

SET.) 

3 

117 

69 

GO  TO  228 

3 

118 

C49 

3 

119 

C50 

THE  NEXT  SEQUENCE  INTEGRATES  THE  CAPACITY 

TO  GIVE  SURFACE  CHARGE 

3 

120 

C51 

AND  INTERFACIAL  TENSION.  J  IS  THE  INDEX  ABOVE  ZERO  CHARGE  POINT. 

3 

121 

70 

DO  140  I*  1, N 

3 

122 

140 

PHI (I)  -  ET All)  -  ETAZ 

3 

123 

141 

DO  72  J* 1 »N 

3 

124 

71 

IF  (PHI! J  ) )  72,  73,  76 

3 

125 

72 

CONTINUE 

3 

126 

73 

CAPZ  «  CAP(J) 

3 

127 

74 

CHGtJ)  -  0. 

3 

128 

75 

GO  TO  78 

3 

129 

76 

CAPZ  *  CAP(J)  -  ( C API J )  -  CAP( J-l) )*PHI <J)/(PHI 1 JJ-PHI ( J-l ) ) 

3 

130 

77 

CHGtJ)  *  0.5*PHI( J )• ( CAP C J)  ♦  CAPZ) 

3 

131 

78 

CHGIJ-1)  *  0.5*PHIIJ-1)*ICAPZ  ♦  CAP! J-l)) 

3 

132 

79 

CHGZ  *  0. 

3 

133 

142 

J2  *  J-2 

3 

133A 

143 

Jl  *  J*l 

3 

1338 

80 

DO  81  I~  1,  J 2 

3 

134 

81 

CHGII)  *  0.5* I  CAP! I )  ♦  CAP! 1*1) )*(PHI (I) 

-  PHI  1 1 *i ) ) 

3 

135 

82 

DO  83  I*  Jl  , N 

3 

136 

83 

CHGII)  *  0.5*1  CAP! I )  *  CAP! 1-1) ) • I  PH  I  1 1 ) 

-  PHI  1 1 -1 ) ) 

3 

137 

84 

K  *  J-l 

3 

138 

85 

K  *  K-l 

3 

139 

86 

CHGIK)  *  CHGIK)  ♦  CHGIK*1) 

3 

140 

87 

IF  IK  -1)  88,88,85 

3 

141 

88 

K  =  J 

3 

142 

89 

K  =*  K  +  1 

3 

143 

90 

CHGIK)  *  CHGIK)  ♦  CHGIK-1) 

3 

144 

91 

IF  I K-N )  89,92,92 

3 

145 

C53 

3 

146 

C54 

WE  HAVE  THE  CHARGE.  WE  INTEGRATE  AGAIN  TO 

GET  INTERFACIAL  TENSION. 

3 

147 

92 

IF  I  PHI I J ) )  93,98,98 

3 

148 

93 

WRITE  OUTPUT  TAPE  6,  94,  PHIIJ),  J 

3 

149 

94 

FORMAT  I/30H  PHIIJ)  IS  NEGATIVE, AND  EQUALS 

,  F5.3,  5H  J  =  12) 

3 

150 

95 

GO  TO  1 

3 

151 
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98  TENS! J)  »  GMAX  -  5.0*PHI I J ) *CHG(J ) 

99  TENS(J-l)  *  GMAX  -  5.0*PH I( J- I ) *CHG( J-I ) 

100  DO  101  I  ■  1«  J2 

101  TENSU)  *  5.0*(CHG(I)+  CHGI I ♦! ) )*(PHI ( I )  -PHI(I*1J) 

102  DO  103  I  -  Jl  ,N 

103  TENSIIJ  ■  5.0*<  CHGt  I )  ♦  CHG( I- 1 ) ) *(PHI ( I )  -  PHl(l-l)) 

104  K  *  J-l 

105  K  *  K-l 

106  TENS(K)  -  TENS! K+l )  -  TENS(K) 

107  IF  (K-l)  108,  108,  105 

108  K  ■  J 

109  K  ■  K*1 

110  TENS(K)  -  TENS(K-l)  -  TENS(K) 

111  IF  (K-N)  109, 112,112 
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MHEN  THERE  IS  NO  AREA  CALCN,  WE  SKIP  THE  2ND  APPROX.  FOR  TENS.  3 

112  IF  (RHO  -  99.999)  167,  120,  115  3 

115  WRITE  OUTPUT  TAPE  6,  116  3 

116  FORMAT  (/  7H  RHO  IS  F8.3,  28H  WHICH  IS  LARGER  THAN  99.999)  3 

323  GO  TO  1  3 

167  IF  ( ABSF (TEST  -  TENS ( 111-  0.01)  113,113,  168  3 

168  NM  *NM+1  3 

169  GO  TO  125  3 

C  3 

113  WRITE  OUTPUT  TAPE  6,  2  3 

402  WRITE  OUTPUT  TAPE  6,  114  3 

114  FORMAT ( //60H  SECONO  APPROX.  BASED  ON  INTEGRATED  INTERFACIAL  T  3 

1ENSI0N  / /66H  RUN  POT  INT . TENS .--ERG/SQ.CM  AREA  C  3 

2AP  POL-RES  /66H  NO.  VOLTS  ASSUMED  INTEGRATED  SO  MM  3 

3  MF/SCCM  0HM-CM2/)  3 

C  3 

170  WRITE  OUTPUT  TAPE  6,  171 , (NUMl ( I ) ,  NUM21I),  ETA(I),  ATENS(I),  3 

1  TENS ( I ) ,  ARE A(  I  I ,  CAP(I),  RES(I),  I  *  1,  N)  3 

171  FORMAT  (IX, 214,  2X,  F7.3,  3X,  2F8.2,  5X,  3F8.3)  3 

172  WRITE  OUTPUT  TAPE  6,  173,  CAPZ,  N,  NM,  J  3 

173  FORMAT  I//42H  THE  CAPACITY  AT  THE  ZERO  CHARGE  POINT  IS  F8.3/12H  TH  3 

1ERE  WERE  I2.25H  DATA  POINTS  IN  THIS  SET./56H  TO  ACHIEVE  A  PRECISIO  3 
2N  OF  0.01  IN  INT. TENSION  REQUIRED  I3,12H  I TERATI ONS. /4H  J*  12)  3 

308  WRITE  OUTPUT  TAPE  6,  35,  DIAM,  M  3 

C60  3 

C61  THE  NEXT  PART  OF  THE  PROGRAM  CALCULATES  THE  CONTRIBUTION  TO  THE  3 
C62  CAPACITY  FROM  THE  DIFFUSE  (GOUY)  AND  COMPACT  (HELMHOLTZ)  DOUBLE  3 

C63  LAYERS, AS  WELL  AS  THE  INTEGRAL  CAPACITIES.  3 

C64  3 

120  DO  178  1*1, N  3 

178  CAPINT  (I)  *  CHG( I )/PHI( I )  3 

174  F J  *  J  3 

175  DO  177  I  *1 , N  3 

176  FI  *  I  3 

177  PSI(I)  =  0.005*( FJ  -  FI)  3 

181  W  *  11.72  *SQRTF(C0NC*(TEMP*273.16)/298.16)  3 

182  V  *  19.46  *298. 16/( TEMP  ♦  273.16)  3 

C 

183  DO  191  I  *1,N  3 

184  F  =  CAPINT(I)*  (PHim-PSIIH)  -  W*SI  NHF  (  V*PSI  ( I))  3 
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179 
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186 
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185  G  »  -CAP  I  NT (I)  -  W»V*COSHF I V*PSI I  I ) )  3 

186  PSI 1  ■  PSI(I)  -I F/G)  3 

187  IF  (A8SFIPSI1  -  PSKIII  -  0.001)  190*190*  188  3 

188  PSI(I)  -  PSI1  3 

189  GO  TO  184  3 

190  PSI  II I  «  PS  11  3 

191  CAP  IG  ( I )  ■  W»SINHF(V*PSHI)  l/PSI  (  I)  3 

193  NN  -1  3 

192  TEST2  -  CAPIGII)  3 

00  2 05  I»1,N  3 

195  CAPIHII)  »  i./l ( i./CAPINTI  I ) )  -  II. /CAPIGII)))  3 

C 

C69  WITH  THIS  IMPROVED  HELMHOLTZ  INTEGRAL  CAPACITY,  BASED  ON  THE  3 

C70  EXPERIMENTAL  RESULTS,  WE  GET  A  BETTER  VALUE  FOR  PSI.  3 

197  F  •  CAPIH(I)«(PHI(I)  -  PSIIII)  -  W»SINHF I V*PSI 1 1 ) )  3 

198  G  ■  -CAPIHII)  -  W»V»COSHFI V*PSI 11)1  3 

199  PSI1  »  PSII  II  -I F/G)  3 

200  IFIA8SFIPSU  -  PSIIII)  -.0001)  203.203,201  3 

201  PSIII)  «  PSII  3 

202  GO  TO  197  3 

203  PSIII)  >  PSII  3 

205  CAPIGII)  -  W*SINHF (V*PSII I ) l/PSI I  I  >  3 

206  IF  I A8SF I CAPIGI 1 )  -  TEST2 )  -  .01)  209,209,207  3 

207  NN  ■  NN  +1  3 

208  GO  TO  192  3 

C71  3 

209  U  «  228. 5*SQRTF( C0NC»298. 16/ 1  TEMP *273. 16 )  )  3 

210  DO  212  I  *  1, N  3 

211  CAPGII)  -  U*COSHF I V*PS 1 1  I )  )  3 

212  CAPHII)  ■  l./l l./CAPI  I)  -1. /CAPGII))  3 

C 

403  WRITE  OUTPUT  TAPE  6,  2  3 

213  WRITE  OUTPUT  TAPE  6,214  3 

214  FORMAT I //61H  RUN  POT.  RAT. POT.  PSI  DIF. CAP  INT.CAP  CHG  3 

1  TENS/65H  NO.  VOLTS  VOLTS  VOLTS  MF/SQCM  MF/SQCM  MC/SQC  3 

2M  ERG/SQCM  /)  3 

215  WRITE  OUTPUT  TAPE  6,  216,  (NUM1II),  NUM2II),  ETA  1 1 ) ,  PHI (I) ,  3 

1  PSIII),  CAP! I ) ,  CAPINT I  I ) ,  CHG 1 1 ) ,  TENS!  I » ,  I  >  1,N)  3 

216  FORMAT  I  IX, 214,  3F8.3,  4F8.2)  3 

217  WRITE  OUTPUT  TAPE  6,218,  ETAZ,  CHGZ.CHGZ,  CAPZ.CAPZ,  CHGZ,  GMAX  3 

218  FORMAT  I/4H  ECM  4X,  3F8.3,  4F8.2)  3 


C 

404  WRITE  OUTPUT  TAPE  6,  2  3 

219  WRITE  OUTPUT  TAPE  6,  220  3 

220  F0RMATI//57H  RUN  RAT. POT  TOTAL  CAPACITY  HELMHOLTZ  DIFF  3 

1USE/62H  NO.  VOLTS  DIFF.  INT.  DIFF.  INT.  DIFF.  I  3 
2NT . /)  3 

221  WRITE  OUTPUT  TAPE  6,  222,  (NUM1II),  NUM2II),  PHI  II) ,  CAPII),  3 

1  CAPINTII),  CAPHII),  CAPIHII),  CAPGII),  CAPIGII),  1  »  l»N)  3 

222  FORMAT  I  IX, 214,  F8.3,  6F8.2)  3 

223  WRITE  OUTPUT  TAPE  6,  224,  N,  NN  3 

224  FORMAT!// 12H  THERE  WERE  12,  25H  DATA  POINTS  IN  THIS  SET./5SH  TO  AC  3 

1HIEVE  A  PRECISION  OF  .01  IN  THE  INTEGRAL  CAPACI TY/39H  OF  THE  HELMH  3 
20LTZ  DOUBLE  LAYER  REQUIREO  I3,12H  ITERATIONS.)  3 

C 

228  WRITE  OUTPUT  TAPE  9,  2  3 
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225  WRITE  OUTPUT  TAPE  9,  21*  ( ETA( I ) ,  CAPII).  NUM(I) ,  I»I,N) 

227  GO  TO  I 

END( l, I, 0,0, 0,0, I, I, 0,0, 0,0, 0,0,0) 


325 

326 
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FORMS  FOR  INPUT  DATA  TO  PROGRAM  3 


The  information  which  must  appear  on  successive  cards,  and  the 
form  in  which  it  must  appear,  is  described  below: 

CARD  1: 

Column  1:  A  "1"  or  "2"  must  appear  in  column  one.  This  signals 
the  beginning  of  a  new  data  set.  A  "l"  sets  the  printer  at  the 
top  of  a  new  page  when  the  output  is  printed.  A  "2"  causes  the 
printer  to  begin  at  the  next  half-page. 

Columns  2-80:  This  space  is  allowed  for  a  title. 

CARD  2: 

Columns  1-8:  RHO,  the  density  of  the  amalgam,  g/cm3  .  Three 
decimal  places. 

Columns  9-16:  HGT,  the  height  of  the  reservoir  above  the  capillary 
tip,  cm.  Three  decimal  places. 

Columns  17-24:  AM,  the  mass  flow  rate,  mg/sec.  Three  decimal 
places. 

Columns  25-32:  DELTA,  the  radius  of  the  cylindrical  or  spherical 
counter  electrode  surrounding  the  drop.  Three  decimal  places. 

Columns  33-40:  TEMP,  the  temperature,  °C.  Two  decimal  places. 
CARD  3: 

Columns  1-8:  GMAX,  the  maximum  interfacial  tension,  erg/cm^ 

(or  dyne -cm).  Two  decimal  places.  This  value  determines 
the  absolute  value  of  the  calculated  interfacial  tension. 

Columns  9-16:  ETAZ,  the  potential  of  zero  charge,  volts.  Three 
decimal  places.  The  potential  is  given  with  respect  to  the  same 
reference  electrode  as  was  used  in  the  capacity  measurements. 
This  value  can  often  be  determined  as  the  zero-current  point 
during  the  capacity  measurements  por  mercury  in 

solutions  with  little  specific  adsorption,  the  zero  charge  potential 
is  -  0.  192  volts  vs.  a  standard  H2  electrode  . 
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Columns  17-24:  CONC,  the  concentration  of  1-1  electrolyte, 
moles/liter.  Five  decimal  places.  For  1-2  or  other  types 
of  electrolytes,  the  last  part  of  the  program  must  be  modified, 
since  the  equations  used  to  calculate  properties  of  the  diffuse 

/Q\ 

double  layer  are  different  ' '  . 

Columns  25-32:  SKAPPA,  the  specific  conductance  of  the  electrolyte, 
ohm  *cm  *  .  Five  decimal  places. 

Columns  33-40:  CRES,  the  resistance  of  the  amalgam  thread  in  the 
capillary,  ohms.  Three  decimal  places. 

CARD  4: 

o 

Columns  1-8:  CAPO,  Kq  in  equation  (6),  fAf/cm  .  Three  decimal 
places. 

Columns  9-16:  CAP1,  in  equation  (6),  pf  cm"  ^  volt"1  .  Three 
decimal  places. 

Columns  17-24:  CAP2,  Kj  in  equation  (6),  plf  cm’^  volt . 

Three  decimal  places. 

This  completes  the  input  of  parameters.  The  next  group  of  cards, from 
l  to  79  in  number,  contains  data  on  drop  lifetime  measured  at  various 
potentials.  The  data  should  appear  in  the  following  form.  The  last  card 
marks  the  end  of  this  data  set. 

DROP  LIFETIME  CARDS: 

Columns  1-8:  ETA1  (M),  potential,  volts.  Three  decimal  places. 
Columns  9-16:  T1  (M),  time,  seconds.  Five  decimal  places. 

LAST  CARD: 

Columns  1-8:  0003.333 

The  next  group  of  cards  contains  the  actual  data.  There  may  be  up  to  79  of 
these.  The  last  card  marks  the  end  of  the  data  set,  and  provides  an  option 
to  terminate  the  calculation  before  performing  the  integration. 
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CAPACITY  DATA  CARDS: 


Columns  1-8:  ETA  (N),  potential,  volts.  Three  decimal  places. 

Columns  9-16:  TIME  (N),  time,  seconds.  Five  decimal  places. 

Columns  17-24:  CAPX  (N),  experimental  value  of  capacity,  micro¬ 
farads.  Three  decimal  places. 

Columns  25-32:  RESX  (N),  experimental  value  of  resistance,  dims. 
Three  decimal  places. 

Columns  33-40:  NUM  (N),  identification  number  consisting  of  up 
to  eight  numerical  digits. 

LAST  CARD: 

Columns  1-8:  0003.333 

This  will  cause  the  calculation  to  run  to  completion,  producing  a  printed 
output  containing  the  capacity,  polarization  resistance,  integrated  values  of 
surface  charge,  interfacial  tension,  integral  capacity,  and  the  capacities  of 
the  Helmholtz  and  diffuse  double  layers. 

ALTERNATE  LAST  CARD: 

Columns  1-8:  0006.666 

This  will  cause  the  calculation  to  terminate  on  completing  the  calculations  of 
capacity  and  polarization  resistance.  In  addition  to  the  printed  output,  cards 
will  be  punched  in  the  correct  format  for  completing  the  integrations  at  a 
later  time. 

In  addition  to  the  direct  path  and  terminated  path  described  above, 
the  program  provides  two  other  alternatives  (Fig.  1). 

MASS  FLOW  INTERPOLATION 

If  the  mass  flow  rate  is  not  constant  throughout  a  series  of  measure¬ 
ments,  a  linear  interpolation  of  mass  flow  rate  can  be  used.  If  this  desired, 
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the  field  (Columns  17-24)  for  AM  on  card  2  should  read  0999.  990.  Card  2 
should  then  be  followed  by  an  additional  card: 

INTERPOLATION  CARD: 

Columns  1-8:  AMI,  initial  value  of  mass  flow  rate,  mg/sec. 

Three  decimal  places. 

Columns  9-16:  AM2,  final  value  of  mass  flow  rate,  mg/sec. 

Three  decimal  places. 

Card  3  should  follow  this  card.  The  rest  of  the  data  input  is  the  same  as 
above.  The  calculation  may  be  terminated  without  integration  by  the 
alternate  last  card. 

INTEGRATION  FROM  PRE-CALCULATED  CAPACITY  DATA: 

If  the  calculation  has  been  terminated  without  integration  in  a  previous 
run,  or  if  data  from  several  calculations,  or  from  the  literature,  are  to  be 
combined  for  integration,  the  first  part  of  the  calculation  can  be  skipped  over. 
If  this  is  desired,  the  data  input  should  be  in  the  following  form: 

CARD  1: 

Column  1:  A  "1"  must  appear  in  column  1. 

Columns  2-80:  Title 
CARD  2: 

Columns  1-8:  This  field  should  read  0099.999. 

Columns  9-32:  These  may  be  blank  or  may  contain  the  same  data 
as  for  the  main  program.  Numbers  appearing  in  these  fields 
will  not  be  used  in  the  calculations. 

Columns  33-40:  TEMP,  the  temperature,  °C.  Two  decimal  places. 
CARD  3: 

9 

Columns  1-8:  GMAX,  the  maximum  interfacial  tension,  erg/cm  . 
Two  decimal  places. 
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Columns  9-16:  ETAZ,  the  potential  of  zero  charge,  volts.  Three 
decimal  places. 

Columns  17-24:  CONC,  the  concentration  of  1-1  electrolyte, 
moles/liter.  Five  decimal  places. 

Columns  25-40:  These  may  be  blank,  or  may  contain  the  same 
data  as  the  main  program.  Numbers  appearing  in  these  fields 
will  not  be  used  in  the  calculations. 

The  fourth  card,  containing  parameters  for  the  trial  capacity  function, 
(equation  6)  is  not  needed,  and  should  be  left  out. 

CAPACITY  DATA  CARDS: 

Columns  1-8:  ETA  (N),  potential,  volts.  Three  decimal  places. 
Columns  9-16:  blank 

Columns  17-24:  CAP  (N),  capacity,  fxf/cm^  .  Two  decimal  places. 
Columns  25-32:  blank 

Columns  33-40:  NUM  (N),  identification  number  consisting  of  up 
to  eight  numerical  digits. 

LAST  CARD: 

Columns  1-8:  0003.333 

With  this  set  of  input  data,  the  program  will  be  entered  at  the  point  where 
the  integrations  are  to  be  performed,  and  will  skip  the  calculations  of  the 
area  and  capacity.  If  it  is  possible  to  obtain  a  complete  set  of  capacity 
data  for  integration  directly  from  the  experimental  results,  this  should  be 
done.  If  the  capacity  is  calculated  in  one  run,  and  the  integrations  are 
performed  in  the  next  run,  some  loss  of  accuracy  results.  This  inaccuracy 
occurs  because  the  drop  area  is  calculated  using  only  the  trial  function  for 
interfacial  tension,  not  the  values  obtained  by  integrating  the  capacity.  In  a 
typical  case,  dividing  the  calculation  this  way  would  introduce  an  error  of 
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about  one  percent  in  the  calculated  capacity  values,  and  a  larger,  cumula¬ 
tive  error  in  the  quantities  obtained  by  integration  which  might  be  as  large 
as  5%. 

Examples  of  input  and  output  data  follow  . 
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RHC  * 

11.660 

GMAX 

’  459.50 

tAPO  « 

41.500 

HGT  * 

45.0C0 

ETAZ 

=  -0.450 

CAP1  » 

114.000 

FLCW  = 

999. 99C 

CONC 

=  O.ICOOO 

"  CAP 2  * 

134.000 

OELTA- 

1.600 

KAPPA 

i  *  0.0385C0 

TEMP  » 

25.00 

CRES  * 

-0. 

THE  MASS 
AMI  ■ 

FLCW  RATE 
12.550  ANC 

USCO  IS 
AM2  ■ 

A  L  INEAR  INTERPOaTION  BETREE7T  " 
11.380 

DRCP  LIFETIME  DATA 

PCT  TIME  GAMMA  DIA 

-1.100  0.27B0C  416.39  0.00258 

-1.000  0. 28 ICC  429.30  0.00250 
-0.900  0.278CO  438.37  0.CC240 

-0.800  0.30000  445.35  0.00252 

-0.700  0.308C0  451.16  0.C0252 

-0.600  O.320CO  455.94  0.00256 

-0.500  O.330C0  459.03  0.00259 

-0.450  0.299C0  459. 5u  0.00231 

AVERAGE  DIAMETER  OF  CAPILLARY  *  0.CC250  MILLIMETERS 

"■NUMBER  Hr  MEASUREMENTS  »  8  - - 
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BASFC  ON  ASSUMED  INTERFACTAL 


■RUN  "  POT  TIME  '  "EXP-CAP  EXP-RE'S  AREA  ^ O^'^m-HTST" 

NO.  VOLTS  SEC  M-FARACS  OEMS  SO  MM  MF/SQCM  0HM-CM2 

770661?  -1.10C  0.21400  0.300  224. 0C0  1.797  16.694 - 3VU6"6- 

7706613  -1.000  0.13400  0.200  160.000  1.293  15.462  1.253 

7706614  -0.90C  0.21700  0.300  BO. 000  1.783  16.874 - C.47I ~ 

7706615  -0.80C  0.20000  0.30C  84.0CC  1.672  17.942  C.478 

7706616  -0.70C  0.15700  0.300  90.000  1.404""Tr*T6Z - 1T.4T3" - 

7706617  -0.600  0.25900  0.600  68.000  1.961  3C.602  C.332 

7706618  -0.50C  0.20300  0.6C0  78.0C0  1.645  36.465  ‘  GV364' 

7706619  -0.45C  0.15700  0.6C0  R2.0CO  1.368  43.869  C.282 
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SECOND  APPROX.  eASEC  CN  INTEGRATED  INTORFAC IAL  TENSION 

RUN  PCT  INT.TENS.~ERG/SQ. CM  AREA  CAP  POL-RES 

NO.  VCLTS  ASSUMED  INTEGRATED  SO  MM  MF/SQCM  0HM-CM2 


7706612 

-1.  lOt. 

A  16 . 39 

AC2.36 

1.798 

16.688 

3.067 

7706613 

-l.CCC 

A29.3Q 

A16.AA 

1.29A 

15.A52 

1.25A 

77066  1A 

-0.90C 

A33.37 

A28.91 

1. 78A~ 

16.820 

“CTA71 

7706615 

-o.eoc 

AA5.35 

A  39 . 71 

1.672 

17.9A1 

C.A78 

7706616 

-0.70C 

A51. 16 

A A  8 . 65 

l.A'OA 

21.363 

C.A13 

7706617 

-0.600 

A55.9A 

A55.31 

1.960 

3C.608 

C.332 

7706618 

-0.560 

A59.03 

A59.0C 

1 . 6A5 

36.‘A72~ 

C.36A 

77066  19 

-0.A5C 

A59.5C 

A59.5C 

1.367 

A3. 878 

C.2R2 

THE  CAPACITY  AT  THE 

ZERO  CHARGE  POINT  IS 

A3. 878 

THERE  WERE 

8  CATA 

PCINTS  IN 

THIS  SET. 

TO  ACMIEVF  A  PRECISION  OF  0.01  IN  INT. TENSION  REQUIRED  2  ITERATIONS. 
J=  8 

AVERAGE  DIAMETER  CF  CAPILLARY  =  0.CC253  MILLIMETERS 

NUMBER  CF  MEASUREMENTS  =  8 


CAPACITY  AMALGAM  8  IN  30. A  AT.IKC  RCOK  77  PAGE  661 


RUN  PCT.  RAT. PCT.  PSI  CIF.CAP  INT. CAP  CFG  TENS- ‘ 


NC. 

VCLTS 

VCLTS 

VCLTS  MF/SQCM 

MF/SCCM 

MC/SCCM 

ERG/SCCM 

7706612 

-1. ICO 

-0.650 

-0.108 

16.69 

22. 9C 

- 1  A  .'89 

A07.36 

7706613 

-l.CCO 

-C .550 

-0.102 

15. A5 

2A.  IA 

-13.28 

Alfc.AA 

77066  14 

-0.9C0 

-0 .A  50 

-0.096 

16.82 

25.92 

-1 1.66 

A  2  8 . 9 1 

T706o  15 

-O.PCC 

-0.360 

-  0 .  C  8  8 

17. 9A 

28.36 

-9.93 

A39.71 

7706616 

-0.7C0 

-0.250 

-0.077 

21.36 

31 ,8A 

-7.96 

A A8.65 

7706617 

-0 .600 

-0.160 

-C.C6C 

3C.61 

35.75 

-5.36 

A55.31 

7706618 

-O.5C0 

-0.0  50 

-0.027 

36.  A  7 

AC.  17 

-2.01 

A59.0C 

77066  19 

-0.A5C 

-0 . 

0. 

A3.8» 

-C. 

C. 

A59.50 

ECM 

-0.A50 

C. 

0. 

A  3 .  P  6 

A3. 88 

0  . 

A  5 9. 50 

32 


CAPACITY  AMALGAP  8  IN  30*4  AT.1KC  BCQK  77  PAGE  661 


RUN - RA77P0T  TCTACtAPACI  TY  HELMHOLTZ  CIFFUSE 

NC.  VCLTS  CIFF.  INT.  DIFF.  INT.  CIFF.  INT. 
7706612  -0.650” "16. 69  ‘  22.90  17.67  27.46  299.07  138.07 


7706613  -0.550  15.45  24.14 

7706614  -0.450  16.82  . 25.92 

7706615  "0.350_  17194  28.36 

7706617  -0.15C  3C.61  35.75 

7706618  -0.050  "36.47  40.17 

7706619  -0.  43.88  -0. 


16.39  29.65  26e.77  129.96 

18. 1C  32.93  238.62  121.77" 

19.65  37.88  206.58  112.87 

-24.41  46.15  ~T7I72T~  101.74— 

40.32  59.50  127.10  89.58 

65.57  86.00  87. W  75.40  ‘ 

111.71  -C.  72.26  C. 


- THERE’ WERE  "8  L'AIA  PL  iNrS  TN  THIS  SET. 

TC  ACHIEVE  A  PRECISION  OF  .01  IN  THE  INTEGRAL  CAPACITY 
‘OF  THE  HELMHOLTZ"  COURT  LAYER  RFQU  TREE  ‘  ~5  TTE  k  A  T  TC  NS.” 
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APPENDIX  II 


Figure  2  -  oscilloscope  trace  showing  points  read. 
Figure  3  -  Format  for  cards  punched  by  graph  reader. 
Listing  of  Program  4. 

Text  describing  input  data  to  Program  4. 

Sample  data. 
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FIG.  2  Oscilloscope  trace  showing  points  read. 
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ormat  for  cards  pun 


non 


DR.  JAMES  N.  BUTLER  -  APRIL  1964  PROGRAM  4 

THIS  PROGRAM  TAKES  DATA.  FROM  THE  GERBER  A-D  CONVERTER  AND  MAKES 
THE  CALCULATIONS  TO  GET  CURRENT  DENSITY  AND  OVERVOLTAGE. 
DIMENSION  KX ( 8  )  ♦  KY (81*  X ( 8 ) •  Y(B>*  U(8>*  T I  ME ( 8 ) *  AREAI8).  ETAI8) 
DIMENSION  CD  (  8  )  •  CDOI8).  T  I  TLF.  (20)  .  NUK(3),  NUMB  (31*  NUM1(3) 

1  RF AD  ini,  DFLTA,  COND 

35  READ  111.  TITLE 

36  TYPE  111.  TITLE 

37  PUNCH  111.  TITLE 

2  READ  102.  RHO 

5  N  =  0 

77  DO  78  I  =  1.3 

78  NUM1 ( I )  =  0 

79  N1  =  0 

6  READ  105.KX(  1  )  .KY(  1  )  .KX  (2)  »KY(2  )  »KX(  3)  *  A.Y  (  3  )  .KX  (4)  »KY ( 4 )  »KX(5)  . 
1KY( 5 ) .KX (6 ) .KY (6  )  »kX ( 7 )  .KY(7 > .KX ( 8 ) »<Y( 8 > *  BIAS.NUM.  M 

38  IF  (M  -  8)  50.  39.  34 

39  KYV1  =  KY ( 1  ) 

40  KYV2  =  KY ( 2 ) 

41  VSC A  =  BIAS 
43  GO  TO  6 

34  IF  IKXI1I  -  998  )  45.  62  .  222 

62  TYPF  108.  NUM1 .  Ml 

63  IF  ( KX (  1  )  -  998  '  1.1.222 

45  KXH 1  =  KX ( 1  ) 

46  KXH2  =  KX ( 2 ) 

47  HSCA  =  BIAS 

49  GO  TO  6 

50  I  F ( M )  11 ,  11  ,  14 

11  IF  ( K  X  (  1  )  -  99P,)  14,  12.  12 

12  TYPc  198,  NUM 1  ,  N1 

13  IF  (  KX ( 1 )  -  993  )  35,  35,  222 

14  N  =  N+l 

7  READ  104,  CURR,  AY,  NUMB.  Mb 

8  IF  ( NUM ( 3 )  -  NUMB (31)  9,  54,  9 

9  TYPE  197,  MUM,  NUMB 

51  GO  TO  6 

54  IF  (V  -  NR)  55,  16,  55 

55  TYPP  197,  M,  MB 
IP  GO  TO  6 

16  YV  =  KY VI  -  KYV2 

17  XH  =  KXH2  -  KXH 1 
82  DO  80  I  =1.  3 
89  NUM1  (I)  =  NIJM  (  I  ) 

81  N 1  =  N 

3  9  PUNCH  199,  NUN' ,  CURR 

18  DO  29  I  =  2,  8 

70  IF  (KX(I))  65.  85,  71 

71  IF  (KY ( I ) )  85.  85,  19 

85  K  =  I  -  1 

86  GO  TO  31 

19  Y (  I )  =  KY ( I  )  -  KYI  1  ) 

29  X(I)  =  KX(I>  -  KXI1I 

21  U(I)  =  BIAS  +  (VSC''  *  Y(T)/  YV) 

22  TIME  (I)  =  HSCA  *  X  (  I  ) /  XH 

23  AREA ( I )  =  4,83597  *  ((AM  *  .001  #  TI'E(I)  /  RHO ) ** ( 2 . / 3 • ) 1 

24  RAD  =  0.28209b  *  SORT  F  ( ARE  A (  1  )  ) 

28  AREA  (  I  )  =  100.  *  AREA!  I  ) 

25  CD( I  1  =  CURR*  .0001/  AREA (  I  ) 
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26  V I R  =  CD ( I )  *  DFLTA  *  RAO  /  ( COND* ( RAD  +  DELTA)) 

27  ETA ( I )  =  U ( I )  -  VIR 

29  CDOHI  =  -0.434294  *  LOGF  {  CD  t  I  )  )  +  8.543*ETA(I) 

87  K  =  8 

31  DO  53  I  =  2.  K 

53  PUNCH  106.  CD(1).  ETA (I).  NUM.  M 
66  DO  32  I  =  2.K 
65  CD< 1 )  =  1000.  *  CDt I ) 

32  PUNCH  110.  U  t I  ) .  TIME!  I).  ARE  A (  I  )  .  ETA ( I ) ♦  CD ( I )  •  CDO(I) 

33  GO  TO  6 
222  CALL  EXIT 

101  FORMAT  C  F  8  •  "4  *  F8.5) 

102  FORMAT  (F8.3) 

104  FORMAT  (F10.4.  F8.4,  51X,  313.  Il> 

105  FORMAT  (13. 13. IX. 13. 13. IX. 13, 13. IX, I  3.13. IX. 13. I  3. IX.  13. 13. IX. 13. I 
13. IX.  13,  I3.5X.F3. 4, lx, 313,  ID 

106  FORMAT  (F10.8,  F8.4,  51X,  313.  ID 

107  FORMAT ( /40H  THE  I. D. NUMBER  ON  THE  PICTURE  CARD  IS  313/  36H  THE  NU 
1MBER  ON  THE  CURRENT  CARD  IS  313/  19H  THEY  DO  NOT  AGREE.) 

108  FORMAT ( /31H  THE  DATA  SET  ENDING  WITH  RUN  313.  13H  IS  FINISHED./ 
114H  IT  CONTAINED  13.15H  PICTURE  CARDS.) 

109  FORMAT ( /24h - RUm  NUMBER  3  I  3 » 9X , F 1 0 . 3 . 10H  MI CROAMP. //48H 

1  U  TIME  AREA  ETA  CURR  DENS  ROUGH  /  50H 

2  VOLTS  SEC.  SO  MM  VULTS  .-iA/SO  C<-»  -LOG  10  /) 

110  FORMAT  (F8.3.F10.5,  F6.3,  F8.3,  F10.4,  FO .  3 ) 

111  FORMAT  (2^44) 

END 
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FORMS  FOR  INPUT  DATA  TO  PROGRAM  4 


The  format  in  which  the  cards  are  punched  is  shown  in  Fig.  3. 

The  vertical  calibration  is  identified  by  an  8  in  column  79;  the  horizontal 
calibration  is  identified  by  a  9  in  column  79;  the  other  columns  of  the 
identification  number  (columns  70  to  78)  are  used  to  identify  the  picture 
from  which  the  data  was  read. 

The  scale  factor  in  the  vertical  calibration  is  the  actual  potential 

difference  in  volts  between  two  calibration  marks,  ytQp  and  yjjQttom  • 

The  scale  factor  in  the  horizontal  calibration  is  the  actual  time  difference 

in  seconds  between  two  calibration  marks,  x,  -  and  x  .  . 

left  right 

The  third  type  of  card  contains  the  actual  data  read  from  the  potential - 
time  curve.  The  first  point  is  a  reference  point,  taken  at  the  bias  voltage, 
at  the  time  corresponding  to  the  initial  point  of  the  curve  (triggering  point 
in  Fig.  2);  the  remaining  7  points  may  be  anywhere  on  the  curve.  The  bias 
voltage  (potential  of  the  base  line)  is  also  included  on  this  card. 

The  remaining  data  required  to  do  the  calculations  is  summarized 
below.  This  is  the  order  in  which  the  cards  must  appear: 

CARD  1; 

Columns  1-8:  DELTA,  the  diameter  of  the  cylindrical  screen 
surrounding  the  drop,  cm.  Three  decimal  places. 

Columns  9-16:  COND,  the  specific  conductance  of  the  electrolyte, 
ohm  cm  ^  .  Five  decimal  places. 

CARD  2: 

Column  1:  A  "1"  must  appear  in  column  1. 

Columns  2-80:  This  space  is  for  a  title. 

CARD  3: 

Columns  1-8:  RHO,  the  density  of  the  amalgam,  g/cm^  .  Three 
decimal  places. 

CARD  4: 

A  vertical  calibration  card  (See  Fig.  3).  "8"  must  appear  in 
column  79. 
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CARD  5: 


A  horizontal  calibration  card  (See  Fig.  3).  "9"  must  appear  in 
column  79. 

CARD  6: 

A  picture-curve  card  (See  Fig.  3).  Column  79  may  contain  any 
digit  from  1  to  7. 

CARD  7: 

Columns  1-10:  CURR,  current,  microamperes.  Four  decimal 
places. 

Columns  11-18:  AM,  mass  flow  rate,  mg/sec.  Four  decimal 
places. 

Columns  70-79:  Run  identification  number.  This  number  must 
agree  exactly  with  the  corresponding  number  on  the  picture 
card  (card  6). 

There  may  follow  any  number  of  cards  of  type  6,  each  followed  by  one  of 
type  7. 

The  horizontal  or  vertical  calibration  is  changed  by  inserting  a  new 
card  of  type  4  or  5  before  the  next  card  of  type  6. 

A  new  set  of  data  is  started  by  a  card  containing  "998"  in  columns 
1-3,  followed  by  cards  of  type  2,  3,  4,  5,  6,  and  7. 

Ordinarily,  the  diameter  of  the  screen  and  the  conductance  of  the 
solution  are  constant  throughout  any  set  of  calculations.  However,  if  it  is 
desired  to  change  the  value  of  DELTA  and  COND,  the  card  before  the  new 
set  of  data  should  contain  "998"  in  columns  1-3  and  "9"  in  column  79. 

This  should  be  followed  by  cards  of  types  1,  2,  3,  4,  5,  6,  and  7. 

The  last  card  of  the  calculation  should  contain  "999"  in  columns  1-3. 
This  program  was  written  for  an  IBM  1620  computer  which  produced 
only  punched-card  output.  It  was  necessary  to  separate  the  output  for 
printing  (PUNCH  statement  32)  from  that  which  was  provided  for  input  to 
Program  5,  which  performs  a  least-squares  analysis  on  the  data  (PUNCH 
statement  53).  This  is  facilitated  by  the  numbers  punched  in  columns  70-79. 
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MAY  13.  1964  PROGRAM  4 
DATA 


1.6  .0385 


DELTA.  COND  FOR  0.1  M  HCL04 


1  THALLIUM  AMALGAM  NO.  8  5-5-64  BOOK  96  PAGES  43  RUNS  1-21.  41.16  MOLE  PCT. 

12.715  RHO  FOR  AMALGAM  8  THALLIUM 


495  213 

063  622 

032284  047310  069260  112213  196175  265156  331143  419130 

130.  3.675 

436285  450314  471260  500225  538198  592176  645161  736143 

131.  3.675 

391286  411294  442240  502199  563175  640158  717146 
118.  3.675 

032290  069298  108262  153240  210221  293201  368189  478172 
60.  3.675 

248  599 

390282  436323  483292  543273  6C3256  695241 
7.02  3.675 

999 


.2000 

.4000 

001.3000 

001.3000 

001.2900 

001.2200 

.5000 

001.0500 


0960430218 
0960430209 
0960430021 
96  43  21 

0960430032 
96  43  32 

0960430  42 
96  43  42 

0960430061 
96  43  61 

0960430199 
0960430092 
96  43  92 
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MAY  13,  1964 
PRINTED  OUTPUT 


PROGRAM  4 


u 

TIME 

AREA 

ETA  CURR  DENS 

ROUGH 

VOLTS 

SEC. 

SO  Mi'-. 

VOLTS  MA/SO  CM 

-LOG  10 

—  RUN  NUMBER 

96  43  2 

130.000 

MICROAMP 

1.318 

.01073 

.102 

1.288 

126.3783 

11.909 

1.282 

.02647 

.187 

1.261 

69.2254 

11.933 

1.249 

.05724 

.314 

1.232 

41.4006 

11.914 

1.222 

.11735 

.506 

1.209 

25.6549 

11.923 

1.209 

.16672 

.640 

1.197 

20.3000 

11.922 

1.200 

.21395 

.756 

1.169 

17.1904 

11.924 

1.190 

.27692 

•  oVB 

1.180 

14.4741 

11.927 

—  RUN  NUMBER 

96  43  3 

131.900 

MICROAMP, 

1.320 

.01001 

.098 

1.290 

133.3448 

11.896 

1.282 

.02504 

.180 

1.259 

72.3906 

11.903 

1.257 

.04579 

.270 

1.239 

46.4106 

11.9U1 

1.238 

.07298 

.369 

1.222 

35.4606 

11.895 

1.222 

.11162 

.490 

1.209 

26.7287 

11.902 

1.212 

.14955 

.595 

1.199 

21.9936 

11.907 

1.199 

.21466 

.757 

1.188 

17.2841 

11.915 

—  RUN  NUMBER 

96  43  4 

118.000 

MICROAMP 

1.295 

.01431 

.124 

1.271 

94.6932 

11.884 

1.257 

.03649 

.232 

1.239 

50.7331 

11.884 

1.228 

.07942 

.390 

1.2)4 

30.2081 

11.896 

1.211 

.12307 

.523 

1.199 

22.5590 

11.893 

1.199 

.17817 

.669 

1.188 

17.6281 

11.909 

1.190 

.23327 

.801 

1.181 

14,7296 

11.922 

—  RUN  NUMBER 

96  43  6 

60.000 

MICROAMP 

■1.225 

.02647 

.187 

1.215 

31.9501 

11.880 

1.200 

.35438 

.303 

1.192 

19.7726 

11.889 

1.184 

.08658 

.413 

1.177 

14.5016 

11.900 

1.171 

.12737 

.535 

1.163 

11.2114 

11.904 

1.156 

.18676 

.690 

1.151 

2.6865 

11.899 

1.148 

.24042 

.817 

1.143 

7.3403 

11.903 

1.136 

.31914 

.967 

1.131 

6.0774 

11.886 

—  RUN  NUMBER 

96  43  9 

7.020 

MICROAMP 

1.079 

.06552 

.343 

1.078 

2.0430 

11.900 

1.067 

.13247 

.540 

1.056 

1.2777 

11.918 

1.043 

.21794 

.765 

1.043 

.9169 

11.948 

1  .0  3  1 

.30341 

.934 

1.031 

.7354 

11.941 

l.u20 

.43447 

1.212 

1.020 

.5788 

11.955 
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MAY  13 

,  1964  PROGRAM  4 

PUNCHED 

OUTPUT 

12637832 

1.2889 

96 

43 

21 

06922540 

1.2611 

96 

43 

21 

04140062 

1.2328 

96 

43 

21 

02565493 

1.2094 

96 

43 

21 

02030C01 

1.1974 

96 

43 

21 

01719040 

1.1892 

96 

43 

21 

01447414 

1.1808 

96 

43 

21 

13334480 

1.2901 

96 

43 

32 

07239068 

1.2598 

96 

43 

32 

04841069 

1.2391 

96 

43 

32 

03548088 

1.2226 

96 

43 

32 

02672872 

1.2091 

96 

43 

32 

02199364 

1.1997 

96 

43 

32 

01728412 

1.1884 

96 

43 

32 

09469320 

1.2713 

96 

43 

42 

05073318 

1.2393 

96 

43 

42 

03020810 

1.2146 

96 

43 

42 

02255900 

1.1994 

96 

43 

42 

01762812 

1.188b 

96 

43 

42 

01472969 

1.1811 

96 

43 

42 

,03195018 

1.2136 

96 

43 

61 

01977268 

1.1922 

96 

43 

61 

,01450166 

1.1777 

96 

43 

61 

,01121144 

1.1651 

96 

43 

61 

,00868655 

1.1516 

96 

43 

61 

,00734032 

1.1435 

96 

43 

61 

,00607740 

1.1319 

96 

43 

61 

,00204304 

1.0782 

96 

43 

92 

,00127779 

1 .0564 

96 

43 

92 

,00091690 

1  .0430 

96 

43 

92 

,00073541 

1.0310 

96 

43 

92 

,00057887 

1 .0204 

96 

4  3 

02 
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APPENDIX  III 


Listing  of  Program  5 . 

Text  describing  input  data  for  Program  5. 
Sample  data. 
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C  JAMES  N  BUTLER.  TYCO  LAOS.  APRIL  1964  PROGRAM  5 

C  HYDROGEN  OVERVOLTAGE,  FRUMKJN  CORRECTION 

C  LEAST  SQUARES  EVALUATION  OF  EXCHANGE  CURRENT  AND  TRANSFER  COEFFICIENT 

DIMENSION  ETA  (100),  CD(IOO),  Xli;"').  Y(100) 

DIMENSION  N1 ( 100) ,  N2 ( 100) ,N3 ( 100) ,  T1TLE120).  NM(IOO) 

1  READ  2,  TITLE 

3  TYPF  2,  TITLF 

4  PUNCH  2.  TITLE 

5  READ  6,  TEMP,  CONC,  ETAZ,  CAPO,  CAP1,  CAP2 
7  PUNCH  8, TEMP,  CONC,  ETAZ,  CAPO,  CAP1,  CAP2 

11  PUNCH  12 
C 

10  X  =  0 

13  X  =  X  +  1 

14  READ  15,  CD ( K  )  ,  ETA ( X ) ,  Nl(<),  N2(X),  N3(K),  NM(K) 

16  IF  ( CD(K ) )  19,  19,  13 

19  K  =  X  -  1 
C 

23  W  =  5039. /(TEMP  +  273.16) 

24  V  =  5892.5/ ( TEMP  +  273.16) 

25  Z  =  0.34C95  *  SQRTF  (CONC  *  (TEMP  +  273.16)) 

C 

26  DO  50  I  =  1,K 

28  PHIO  =  -  F  T  A (  I  )  +  ETAZ 

2  9  CAP  =  CAPO  +  ( CAP  1  *  Pill  0 )  +  (CaP2  *  PH  10*  PHIO) 

30  P  =  PHIO  *  CAP 

C23  THE  PSI-EQUAIION  IS  SOLVED  bY  1TERAT1U.,  USING  NEWTON'S  METHOD 

31  PSI  =  -0.05 

32  F  r  CAP  *  PSI  -  P  +  7  * ( EXPF  (V*PSI)  -  EXPF  (-V*PSI>) 

33  G  =  CAP  +■  Z*V* ( FXPF  (V*PrI)  +  EXPF  (-V*PSI)) 

34  H  =  PSI  -  ( F/G) 

35  I F  (  ( AtiSF  ( H  -  PSI))  -  .0001  )  37,37,36 

36  PSI  =  H 
GO  TO  32 

C24  FND  OF  ITERATION 

37  PSI  =  H 

27  CDO  =  -  0,4  342°4  *  LOGF(CD(I)>  +  0.5  *  W  *  ( ET  A (  I  )  -  PSI) 

38  X ( I )  =  ETA ( I )  +  PSI 

39  Y  (  I  )  =  PSI  +  0.434294*  (  LOGF  (CD(  I  )  )  ) /.. 

21  CD ( I  )  =  1009.  «  CP(  I  ) 

4  PUNCH  41,  N1(I),  N2  (  I  )  ,  N3  (  I  >  »  EIA(I),  C  D  (  1  )  ,  CAP,  PSI,  PHIO,  CDO, 

IX ( I  ) ,  Y  (  I  ) 

50  CONTINUE 

C 

51  PUNCH  ■=.?,  X 
C 

C 26  NOW  COMES  THC  LEAST-SQUARES  PART 
53  SX  =  0. 

SY  =  0. 

SXX  =  0. 

SX Y  =  0. 

SY Y  =  0. 


5  4 

II 

•— 1 

c 

vO 

c 

Cl 

1,  X 

55 

SX  =  sx  + 

X  (  I  ) 

56 

SY  =  SY  + 

Y  (  I  ) 

57 

SXX  =  SXX 

+  X (  I  )*X (  I  ) 

58 

SX Y  =  SXY 

+  X (  I )*Y (  I  > 

59 

SY Y  =  SY  Y 

+  Y(  I  )*Y(  I  ) 

60 

CONT I NUE 
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61  FN  =  < 

62  DN  =  l./FN 

63  AVX  =  SX  *  DN 

64  AVY  =  SY  *  DN 

65  DNN  =  1 • / ( FN  -  l.» 

66  S  =  SXX  -  SX  *  AVX 

67  COVX  =  S  *  DNN 

68  COVY  =  (SYY  -  SY  *  AVY)  «  DNN 
6<?  COVXY  =  ( SXY  -  SX  *  AVY)  *  DNN 

70  STDX  =  SQRTF  (COVX) 

71  STDY  =  SQRTF  (COVY) 

72  CORR  =  COVXY/ (STDX  *  STDY) 

73  A  =  COVX Y/COVX 

74  B  =  AVY  -  A#AVX 

75  SR  =  SQRTF((FN-1.)*(C0VY-A»A»C0VX)/(FN  -  2.)) 

76  SI  =  SR  *  SQRTF ( DN  +  AVX  #  AVX  /  S) 

77  SS  =  SR  *  SORTF  (l./S> 

78  PUNCH  79  *  CORR.  A.  B*  SR*  SS.  SI 

80  CDOP  =  -B  *  V 

81  CONFA  =  1.65  *  SS 

82  CONFB  =  -1.65  *  SI  »  W 

83  PUNCH  84.  A,  CONFA.  CDOP.  CONFB 

C27  THUS  ENDS  THC  LEAST  SQUARES  PART 

85  READ  2,  TITLE 

86  IF  (TITLE(l)  -  .455544)  87,  88,  87 

88  TYPF  2,  TITLF 

89  PUNCH  2,  TITLF 

90  CALL  EXIT 

87  GO  TO  3 
C 

2  FORMAT  ( 20A4 ) 

6  FORMAT  (F8.2,  F8.5,  4F8.3) 

8  FORMAT  ( / /8H  TEMP  =  F8.2/8H  CONC  =  F8.5/8H  ETAZ  =  F8.3/8H  CAPO  = 
1F8.3/8H  CAP  1  =  F8.3/8H  CAP2  =  F8.3//37H  LF AST-SOUARFS  FOR  ALPHA  AN 
2D  I O-PRIME  ) 

12  FORMAT  ( 7 7H 1  RUN  NO.  ETA  CURR  DENS  CAP  PSI  PHI  - 

1  LOG  I  OPR  I  HE  X  Y  /62H  VOLTS  MA/SQ.CM  MF/SQ. 

2CM  VOLTS  VOLTS  (ALPHA=.5>  ) 

15  FORMAT  (F10.8,  F8.4,  51X,  313,  ID 
41  FORMAT  (IX,  313,  F8.3,  F10.4,  3F8.3,  3F9.3) 

52  FORMAT  (/12H  THERE  ARE  I2.24H  DATA  POINTS  IN  THIS  SET) 

79  FORMAT  (//31H  CORRELATION  BETWEEN  X  AND  Y  =  F10.5/22H  SLOPE  OF  LIN 

IE  Y ( X )  =  F10.5/26H  INTERCEPT  OF  LINE  Y l X )  =  F10.5/31H  STANDARD  ERR 
20R  OF  EST I  MAT  t  =  F10.5/31H  CONFIDENCE  INTERVAL  PARAMETERS/22X, 10H 

3  SLOPE  F10.5/22X, 10H  INTERCEPT  F10.5) 

84  FORMAT  ( / 4 1 H  CATHODIC  TRANSFER  COEFFICIENT,  ALPHA  =  F10.5/32H  95 

1PERCENT  CONFIDENCE  LIMITS  =  F8.3//41H  CORRECTED  EXCHANGE  CURRENT, 
2-LOG  10  =  F9.3/32H  95  PERCENT  CONFIDENCE  LIMITS  =  F8.3) 

END 
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FORMS  FOR  INPUT  DATA  TO  PROGRAM  5 


The  input  to  Program  5,  the  second  part  of  the  overvoltage  calcula¬ 
tions,  should  consist  of  the  following  cards: 

CARD  1; 

Column  1:  A  "1"  should  appear  in  column  1. 

Columns  2-80:  This  space  is  for  a  title. 

CARD  2: 

Columns  1-8:  TEMP,  the  temperature,  °C.  Two  decimal  places. 

Columns  9-16:  CONC,  the  electrolyte  concentration,  in  moles/liter. 
Five  decimal  places. 

Columns  17-24:  ETAZ,  the  zero  charge  point  of  the  amalgam  with 
respect  to  the  same  reference  electrode  as  the  overpotential 
measurements.  Three  decimal  places. 

Columns  25-32:  CAPO,  Kq  in  equation  (6), 
decimal  places. 

Columns  33-40:  CAPl,  in  equation  (6), 

Three  decimal  places. 

Columns  41-48:  CAP2,  ^  in  equation  (6), 

Three  decimal  places. 

This  is  followed  by  data  cards  in  the  following  format: 

DATA  CARDS: 

Columns  1-10:  CD  (I),  current  density,  amp/cm2.  Eight  decimal 
places. 

Columns  11-18:  ETA  (I),  overpotential,  volts.  Four  decimal  places. 

Columns  70-79:  Identification  number 

Program  4  produces  output  cards  in  precisely  this  format,  so  that  it  is 
necessary  only  to  select  the  desired  cards  for  input  to  Program  5. 


|lf/cm2  .  Three 
^lf  cm'2  volt"1  . 
/tlf  cm’2  volt’2  . 
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LAST  CARD: 


The  last  card  of  the  data  set  should  be  blank. 

The  last  card  of  the  entire  deck  should  read  "END"  in  columns  1-3. 


NOTE: 

The  program  has  been  written  with  the  convention  that  cathodic 
overpotentials  are  taken  to  be  positive.  A  zero- charge  point  which  is  nega¬ 
tive  with  respect  to  the  zero-overpotential  point  should  be  given  a  positive 
value  also.  If  both  the  overpotentials  and  the  zero-charge  point  are  given 
the  conventional  signs  for  electrode  potentials  (negative  for  cathodic,  positive 
for  anodic),  the  program  must  be  modified  as  follows: 

Statement  28  should  read  "PHIO  =  ETA  (1)  -  ETAZ". 

A  statement  should  be  inserted  after  statement  37,  which  reads 
"ETA  (1)  =  -  ETA  (1)  ". 

A  statement  should  be  inserted  after  statement  21,  which  reads 
"ETA  (1)  =  -  ETA  (1)"  . 
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1  MERCURY  1-22-64  BOOK  06,  P4f^  10,  RUNS  1-11 

26.2  .1  .165  26.6  16.66  9.66  96  19  MERCURY 


.10711318 

1.2802 

96 

19 

11 

.06614743 

1.2555 

96 

19 

11 

.04126667 

1.2308 

96 

19 

11 

.03089814 

1.2166 

96 

19 

11 

.02365379 

1.2027 

96 

19 

11 

.11971218 

1.1938 

96 

19 

11 

.08170679 

1.2636 

96 

19 

22 

.04558101 

1.2325 

96 

19 

22 

.03012561 

1.2104 

96 

19 

22 

.01888411 

1.1866 

96 

19 

22 

.01380436 

1.1718 

96 

19 

22 

.00924673 

1.1528 

96 

19 

22 

.00696264 

1.1405 

96 

19 

22 

.02097022 

1.1915 

96 

19 

33 

.01066125 

1.1576 

96 

19 

33 

.00635090 

1.1324 

96 

19 

33 

.00385461 

1.1084 

96 

19 

33 

.00290654 

1.0943 

96 

19 

33 

.00255973 

1.0886 

96 

19 

33 

.00476340 

1.1191 

96 

19 

44 

.00270185 

1.0935 

96 

19 

44 

.00164911 

1.0713 

96 

19 

44 

,00220697 

1.0839 

96 

19 

65 

,00079137 

1.0372 

96 

19 

65 

.00046663 

1.0163 

96 

19 

65 

.00036849 

1.0062 

96 

19 

65 

.00024047 

.9824 

96 

19 

76 

.00014127 

.9585 

96 

19 

76 

.10011258 

.9491 

96 

19 

76 

END  OF  RUN 

96  19  MERCURY 
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1  MERCURY  1-22-64  ROOK  96,  PAGE  19,  RUNS  1-11 


TEMP 

= 

26.20 

CONC 

= 

.10000 

ETAZ 

= 

.165 

CAPO 

= 

26.600 

CAP  1 

= 

18.650 

CAP2 

r 

9.650 

LEAST-SQUARES  FOR  ALPHA  AND  IO-PRJME 


1  RUN  NO 

• 

ETA 

VOLTS 

CURR  DENS 
MA/SU.CH 

CAP  PSI 

MF /SQ.CM  VOLTS 

PHI 

VOLTS 

-LOG IOPR I  ME 
( ALPHA=.5 ) 

X 

Y 

96  19 

1 

1.280 

107.1131 

17.802 

-.  1 16 

-1.115 

12.729 

1.163 

-.174 

96  19 

1 

1.255 

66.1474 

17.737 

-.115 

-1.090 

12.719 

1.139 

-.185 

96  19 

1 

1.230 

41.2666 

17.684 

-.114 

-1.065 

12.704 

1.116 

-.196 

96  19 

1 

1.216 

30.8981 

17.659 

-.113 

-1.051 

12.704 

1.103 

-.203 

96  19 

1 

1.202 

23.6537 

17.638 

-.112 

-1.037 

12.696 

1.090 

-.209 

96  19 

1 

1.193 

19.7121 

17.626 

-.112 

-1.028 

12.696 

1.081 

-.213 

96  19 

2 

1.263 

81.7067 

17.757 

-.116 

-1.098 

12.699 

1.147 

-.180 

96  19 

2 

1.232 

45.5810 

17.687 

-.114 

-1.067 

12.676 

1.118 

-.194 

96  19 

2 

1.210 

30.1256 

17.649 

-.113 

-1.045 

12.660 

1.097 

-.203 

96  19 

2 

1.186 

18.8841 

17.618 

-.111 

-1.021 

12.650 

1.074 

-.214 

96  19 

2 

1.171 

13.8043 

1  7.604 

-.110 

-1.006 

12.656 

1.060 

-.221 

96  19 

2 

1.152 

9.2467 

17.593 

-.109 

-.987 

12.661 

1.042 

-.230 

96  19 

2 

1.140 

6.9626 

17.589 

-.109 

-.975 

12.675 

1.031 

-.237 

96  19 

3 

1.191 

20.9702 

17.624 

-.112 

-1.026 

12.649 

1.079 

-.211 

96  19 

3 

1 .  1d7 

10.6612 

17.595 

-.110 

-.992 

12.642 

1.047 

-.227 

96  19 

3 

1.132 

6.3509 

17.589 

-.  108 

-.967 

12.643 

1.023 

-.239 

96  19 

3 

1.108 

3.8645 

17.594 

-.  107 

-.943 

12.647 

1.000 

-.250 

96  19 

3 

1.094 

2.9065 

17.602 

-.106 

-.929 

12.644 

.987 

-.257 

96  19 

3 

1.088 

2.5597 

17.606 

-.  106 

-.923 

12.649 

.982 

-.260 

96  19 

4 

1.119 

4,7634 

17.590 

-.108 

-.954 

12.650 

1.011 

-.246 

96  19 

4 

1.093 

2.7018 

17.602 

-.106 

-.928 

12.669 

.986 

-.259 

96  19 

4 

1.071 

1.6491 

17.623 

-.105 

-.906 

12.686 

.965 

-.270 

96  19 

6 

1.083 

2.2069 

17.610 

-.106 

-.918 

12.671 

.977 

-.263 

96  19 

6 

1.037 

.7913 

17.674 

-.103 

-.872 

12.702 

.933 

-.287 

96  19 

6 

1.016 

.4666 

17.716 

-.  102 

-.851 

12.745 

.913 

-.300 

96  19 

6 

1.006 

.3684 

17.740 

-.  101 

-.841 

12.758 

.904 

-.305 

96  19 

7 

.982 

.2404 

17.803 

-.100 

-.817 

12.732 

.881 

-.315 

96  19 

7 

.958 

.1412 

17.877 

-.099 

-.793 

12.750 

.859 

-.327 

96  19  7  .949  .1125  17.909  -.098 

THERE  ARE  29  DATA  POINTS  IN  THIS  SET 

CORRELATION  BETWFEN  X  AND  Y  =  .99884 

SLOPE  OF  LINE  Y ( X )  =  .50967 

INTERCEPT  OF  LINE'Y(X)  =  -.76365 

STANDARD  ERROR  OF  ESTIMATE  =  .00219 

CONFIDENCE  I NTFRVAL  PARAMF TFRS 

SLOPE  .00472 

INTERCEPT  .00485 

CATHODIC  TRANSFER  COEFFICIENT,  ALPHA  = 

95  PERCENT  CONFIDENCE  LIMITS  =  .007 

CORRECTED  EXCHANGE  CURRENT,  -LOG  10  =  12 

95  PERCENT  CONFIDENCE  LIMITS  =  -.134 

-.784 

50967 

.854 

12.765 

.850 

-.333 
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